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Abstract. 

I discuss the systematic modifications to the perturbative QCD factorization 
approach in high energy i+A, p+A and A+A reactions. These include transverse 
momentum diffusion manifest in the Cronin effect and a small increase in the dijet 
acoplanarity; nuclear size enhanced power corrections that lead to shadowing 
in deeply inelastic scattering and suppression of single and double inclusive 
hadron production at forward rapidity at RHIC but disappear as a function of 
the transverse momentum; inelastic attenuation of the jet cross sections or jet 
quenching that persists to much higher px- 
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1. Introduction 

The interpretation ^ |2] of the copious experimental results from the past, present 
and future heavy ion programs necessitates a reliable theoretical framework for 
the calculation of a variety of moderate and large transverse momentum processes, 
Q 2 ~ Pt ^ Aq CD , that can be systematically extended to incorporate corrections 
arising from the many body nuclear dynamics. Such theoretical framework is the 
perturbative QCD factorization approach py, a direct generalization of the naive 
parton model 0J. According to the factorization theorem, the observable hadronic 
cross sections can be expressed as 

Chadronic = °partonic Zj) pr\ fJ>fi, Pdj) ® | Y\_ fii/hi Pf i) | 

Dhj/jiZj'Pdj)?;, (1) 



n 
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where <8 denotes the standard convolution over the internal kinematic variables of 
the reaction. In Eq. JQ) /i r , [ifi and \idj are the renormalization, factorization and 
fragmentation scales, respectively. <j>i/ht( x h M/«) i s the distribution function (PDF) 
of parton "i" in the hadron hi and D^./j^Zj, fidj) is the fragmentation or decay 
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function (FF) of parton "j" into hadron hj. Factorization not only separates the 
short- and long-distance QCD dynamics but implies universality of the PDFs and 
FFs and infrared safety of the hard scattering partonic cross sections. If the kxi 
and/or the jxi dependence of the PDFs, FFs and <x par tonic is kept explicit, one arrives 
at a kr and/or jx factorized form (BJ- At present, limited statistics hampers any 
reliable extraction of such transverse degrees of freedom. Integrating kxi and jxi out 
leads to the most commonly employed double collinear limit of QCD. 

Factorization has been discussed in detail for a limited number of processes. 
Explicit examples of some fundamental cross sections are given below to lowest order 
and leading twist. All scales are suppressed and the definitions of the kinematic 
variables can be found in [6]. These examples include: 

(i) Electron-positron annihilation into hadrons [7]: 

= N ^ Y,Q* (1 + cos 2 M D hj/j ( Zj ) . (2) 

Here x — 2p ■ q/q 2 , y = ^(1 — cos# cm ) and Qi is the fractional electric charge 
of (anti) quarks in units of e. For a discussion on the structure of the near-side 
correlation function from fragmentation in e + + e~ annihilation see 8 . 

(ii) Deeply inelastic lepton-hadron scattering (DIS) with the longitudinal and 
transverse structure functions: 

^r = i]TQ 2 0/(x) , F L =0. (3) 

i 

In Eq. © x = x B — -q 2 /2p ■ q = Q 2 /2m N v with v = E — E', Recently 
the nuclear enhanced high twist corrections to the DIS cross sections have been 
evaluated within the factorization approach [5] I1U) . Similar results from final 
state multiple scattering have been found in 

(iii) The Drcll-Yan process |12| : 



dvi+i- v^, , u , , 47ra 2 m Q 2 



Xa)4>l/N> M OA e T m T 6 (l 2 - ( X "P + ^p'f) ■ (4) 



In Eqs. J5J and (0J N c = 3 is the number of colors. The momentum fractions 

X a = Pt I Pa and X b = Pb/ P b + - 

(iv) Inclusive hadron production in N + N collisions to leading power and leading 
power corrections JS| • The single and double inclusive distributions accessible to 
0(a 2 ) at leading twist read [T3] : 



da NN _ f dzi , . f 4> a /N(x a ) 



dyid 2 p> 



Ti 



abed 



i 



x a S + U/zi 



x 4 \M ab ^ d \ 2 , (5) 



S Xb 

d^NN _ 8(Atp — 7r) ^ f dzi <Pa/ N (x a ) 



dy 1 dy 2 d 2 p Tl d 2 pT 2 PTiPt 2 Z\ 
a 2 s fo/N&b 



D hi/c(zi) D h2 / d (z 2 ) ■ 



Mab^cdl 2 ■ (6) 
O Xb 

In Eqs. ||SJ| and © S, T and U are the hadronic Mandelstam variables |14|. 
The dependence on the small momentum fraction Xb, relevant to forward rapidity 
studies at RHIC, is explicitly isolated. 
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Figure 1. Left panel: the differential double inclusive tt° distributions versus 
the fragmentation momentum fractions zi, Z2 for a trigger pion of px = 5 GeV 
and associated pion of pj< = 2 GeV. Right panel: the same cross sections for a 
trigger pion of pt = 5 GeV and associated pion of pt = 8 GeV. Bottom panel: 
the anti-correlation between Ztrig and pt assoc f° r a fixed pt trig = 5 GeV. The 
dashed line simulates a trigger bias Az tr ; g = 0.15. 



(v) Heavy quark and heavy quark bound state production |15|. 
1.1. The fragmentation seesaw analogy 

One of the exciting theoretical questions, prompted by the recent experimental 
results ^3 E| j is about the nature of the hadron production mechanism at small and 
moderate pt in A + A reactions |18U19ll!?Ull21| . The measurement of correlations [221 
in p + p and especially d + Au collisions already provide sufficient evidence that back- 
to-back hard partonic scattering and fragmentation, Eqs. JSJ and ©, give a dominant 
contribution to the particle cross sections. Corrections in the basic perturbative 
formulas arise both from higher orders in a s and from higher twist. These may be 
enhanced by the nuclear size and the density of the quark-gluon plasma (QGP) in the 
deconfined phase. Such corrections can be systematically organized in the framework 
of the factorization approach as follows |24| : 

Madron = ^ (l + Cf } a s + Cf + Cf ) a\ + • ■ •) 

(4) 

+ ^-(i + cf)«. + cf^ + cfa» + -)rW 



4 



+ ^ (l + C[^a s + Cf^al + C^al +■■■) T«0 + - - • (7) 

In Eq. © rW are the twist "i" non-perturbative matrix elements. Odd- twist matrix 
elements play an important role in spin physics |25| but are here neglected for 
simplicity. 

As a test of the dominance of the lowest order, leading twist term one can 
experimentally look for a seesaw fragmentation analogy. The seesaw mechanism |26| 
extends the Standard Model |23 by including a massive right-handed Majorana 
neutrino. The physical neutrino masses are given by M„ = MtM^M^. When 
Mr is very large and hence vr is unobservable, M„ can become very small. 

Similar anti-correlation is present in Eq. ©. Naively, one might expect that a 
high-pT hadron trigger in the near side will fix the fragmentation momentum fraction 
Ztrig- In contrast, 6(j>ti/zi — PT2/Z2) implies that ztrig is inversely proportional to 
Pt assoc- Figure n shows the differential distributions da hlh2 / dy\dy2dpT\dpT2dzi for 
different choices of pr trig and pt assoc ■ Experimental verification of the anti-correlation 
presented in the bottom panel of Fig. ^ can provide critical additional evidence in 
support of the dominance of the perturbative QCD hadron production mechanism 
and the cross section hierarchy summarized by Eq. (JJJ. 



2. Transverse momentum diffusion 



Hard processes in the nuclear environment involve multiple parton interactions before 
or after the large Q 2 collision that are sensitive to the properties of the nuclear matter. 
In this section we consider the multiple elastic and incoherent scattering of energetic 
quarks and gluons that penetrate hot and cold QCD matter. Traverse momentum 
diffusion reflects the (p^)-kick per unit length in the medium, given by the transport 
coefficient fi 2 /\ qtg , and is the strong interaction dynamics equivalent of the Moliere 
multiple scattering in QED. 

Approximate solutions for the parton broadening can be obtained to leading power 
and leading power corrections in the large "+" lightcone momentum p + |28| . For the 
instructive case of a normalized forward monochromatic beam (p = pu = P) the 
initial quark or gluon distribution reads: 

d 3 N i 

In Eq. (JSJ the constraint on the parton's "— " lightcone component from the p 2 = 
on-shell condition is also shown. In the small angle approximation, rescattering leads 
to an approximately Gaussian form |28| : 



d 3 N f (p+,p±) 1 e 

dp+d 2 p ± 'p-=^ ~ 2^ x^H 



1 V2 2P 



(9) 



where \ — L/X is the opacity or mean number of scatterings and £ = 0(1). It should 
be realized that this form is not expected to describe well the large angle scattering. 
The tails of the large distributions are expected to be power-law like and die out 
at a much slower rate than the Gaussian in Eq. 10. This may lead to uncertainties 
in the separation of the jet and the underlying event |22) . 
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Figure 2. Left panel from |31l : calculation of the Cronin effect for 7r+ + 7r~ 
at y/sjvjv = 27.4 and 38.8 GeV at midrapidity. The nuclear modification is 
identified from the ratio of the cross sections in p + W and p + Be reactions. Right 
panel from |31| : predicted magnitude and centrality dependence of R^\(y>t) 
at ^/sjviv = 200 GeV and midrapidity at RHIC in the absence of significant 
antishadowing. 



2.1. Applications of transverse momentum diffusion 

From Eq. © it is easy to demonstrate that the the accumulated acoplanarity 
momentum squared by an energetic parton probes a line integral through the color 
charge density: 

V 2 o<» f j SCfl™? 



2£^P±pS(L), static 

2 (10) 

"» 2 Aj_ r ' 

In Eq. (|1(J|) the factor 2 comes from 2D diffusion, and p s is the effective gluon density. 
For the 1+1D Bjorken expansion scenario A± is the transverse area of the interaction 
region, tq is the initial equilibration time and dN 9 /dy is the effective gluon rapidity 
density. 

Dynamical nuclear-induced modification in multi-particle production can be 
studied through the ratio |14j : 



R (n) _ d(jl AB hn I ' ' ' d ynd 2 p Tl [ ' ' d 2 p Trl 

AB ~ (N%°g) da^ h " /d Vl ■ ■ ■ dy n d^p Tl ■ ■ ■ d^p Tn ' 



Centrality dependence is implicit in Eq. (|f I (I . The first and possibly most sensitive 
indicator of multi-parton dynamics in cold nuclear matter is the Cronin effect |29| in 
the single inclusive hadron production in p + A reactions, ^^(pt). 

Elastic transverse momentum diffusion at midrapidity is manifest in the Cronin 
enhancement at pr > I — 2 GeV and suppression for pt < f — 2 GeV. For a survey 
of theoretical models see j^O]. Calculations of the nuclear modification in low energy 
p + A reactions employ /i 2 /A = O.f — 0.f4 GeV 2 and are shown in the left hand side 
of Fig. [21 The predicted |3J centrality dependence of the Cronin effect in d + Au 
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Figure 3. Left panel from I28I : the acoplanarity momentum projection per 
parton ^ |) and the corresponding width <rp ar of the away-side correlation 
function for p + p and d + Au. The approximations of «trig 1 an d constant 



(\ k Ty\) f 



0.75 had been used. Right panel from |22]: preliminary PHENIX 



data on the vacuum and medium induced acoplanarity momentum projection 
per parton in cold nuclear matter for realistic finite z tI i g . Note the absence of 
significant difference between the p + p and d + Au measurements. 



reactions at y = and RHIC ^/snn = 200 GeV is given in the right hand side of 
Fig. [3 and confirmed by experimental data |221 1231 132] . It is important to emphasize 
that even a small broadening may lead to a noticeable Cronin effect since it is amplified 
by the steepness of the partonic/hadronic spectra. 

It is interesting to look for other manifestations of transverse momentum diffusion 
such as increased dijet acoplanarity, which will be reflected in the increased width of 
the away-side correlation function 

c 2 {^)- 1 diVdlJot 



2 , ^Far 



^iNcar /IF 



In Eq. H12|) the near-side width UNoar of C^Ai^s) is determined by jet fragmentation 
and CTFar reflects (ky) to t = (^t)pp + ( k r)nud.- 

The left panel of Fig. shows the upper limit of the medium induced acoplanarity 
relative to the vacuum p+p result. The predicted increase in the width of the away-side 
correlation function is small even in the z tr i g — > 1 limit. For realistic fragmentation 
momentum fractions, see Fig.^ the experimentally observed enhancement of (Ak^) to t 
should be significantly smaller. In addition, for dihadron angular correlations there 
are no amplification effects from the steepness of the underlying parton spectra. The 
right panel of Fig. [31 indeed shows the lack of significant difference between p + p and 
d + Au measurements and certainly excludes monojet-based models. 
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Figure 4. Left panel from all-twist resummed F2[A)/F$(D) calculation 
from versus CERN-NA37 and FNAL-E665 data on DIS on nuclei. The band 
corresponds to the choice £ 2 = 0.09 — 0.12 GeV 2 . Data-Theory, where Ajj_y is 
computed for the set presented by circles, also shows comparison to the EKS98 
scale-dependent shadowing parametrization. Right panel from 9 : CERN-NA37 
data on Fi^Sn) / F2(C) show evidence for a power-law in 1/Q 2 behavior consistent 
with the all-twist resummed calculation. The bottom right insert illustrates the 
role of higher twist contributions to Fl on the example of R = ax / a T ■ 



3. Nuclear enhanced power corrections 

A class of corrections that can be naturally incorporated in the perturbative QCD 
factorization approach is associated with the power suppressed ~ \/Q n , n > 1 
contributions. The higher twist terms in Eq. Q are typically neglected in reactions 
with "elementary" nucleons for Q 2 > 1 GeV 2 . However, in the presence of nuclear 
matter such corrections can be enhanced by the large nuclear size ~ A 1 / 3 . These lead 
to dynamical nuclear shadowing. We discuss this mechanism in detail below. 

Hard scattering in nuclear collisions requires one large momentum transfer Q ~ 
xP 3> h-QCD with parton momentum fraction x and beam momentum P. A simple 
example is the lepton-nucleus deeply inelastic scattering (DIS). The transverse area 
probed by the virtual meson ~ 1/Q 2 is very small for moderate and large Q 2 processes. 
Hence, the struck quark will propagate in the nucleus and interact with nucleons along 
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the same impact parameter. If Q 2 <C m 2 N the transverse size of the wave packet will be 
comparable or larger than the size of the nucleon. In this case the rapid fall-off of the 
parton wave function effectively limits Q 2 > Q\. The effective longitudinal interaction 
length probed by the virtual meson of momentum is characterized by 1/xP. If the 
momentum fraction of an active initial-state parton x <C x c — l/2mjv?"o ~ 0.1 with 
nucleon mass tun and radius ro, it could cover several Lorentz contracted nucleons of 
longitudinal size ~ 2ro(jn,N / P) in a large nucleus. 

Each of the multiple soft interactions is characterized by a scale of power 
correction per nucleon 

?«^^b). (is) 

with the matrix element (p\F 2 \p) = g lim x _jo xG(x, Q 2 ). A large nucleus will provide 
~ A 1 / 3, enhancement of the high twist shadowing contribution. We note the Eq. (f 1 31) 
has the geometric average for minimum bias reactions incorporated in the definition 
of£ 2 . 



3.1. Applications of dynamical power corrections 

In Ref. [5] we resummed the nuclear enhanced high twist corrections and identified the 
modification to the leading twist and lowest non- vanishing order in a s contribution to 
the longitudinal and transverse structure functions: 



T 

F£(x, Q 2 ) » A ff T) (x, Q 2 ) + ^ F$(x, Q 1 i . ( I i 



f T) M 2 ) + f 2 



These results, when compared to Eq. |J3J|, indicate that the essential difference from 
the final state coherent multiple parton scattering is the generation of dynamical mass 
and the consequent rescaling in the value of Bjorken-ir. In addition, there is a novel 
contribution to F£. 

Calculations of nuclear shadowing versus xb with £ 2 = 0.09 — 0.12 GeV 2 are given 
in left panel of Fig. 01 The right panel shows the Q 2 dependence of power corrections 
and the enhancement in the longitudinal structure function, see Eq. (|15f) . reflected in 
R = a L /o- T - 

This resummation, derived within the framework of the pQCD factorization 
approach has definite advantages: 

• It relies on standard PDFs 33; and its particularly easy to implement numerically. 
The final state soft scattering of the struck parton in the medium generates 
dynamical parton mass w^ yn = ^(A 1 ^ 3 — 1) via the coupling to the background 
chromo- magnetic field given by the two gluon correlation function, Eq. (|13fl . The 
necessary additional energy leads to the rescaling of the value of Bjorken-ir. This 
physical interpretation becomes transparent when we examine the behavior of 
v + A DIS with charged current exchange and charm quark in the final state 10 
with large physical mass M c . Such processes are allowed via the CKM matrix 



mixing and we find: 



M 2 ^(A 1 / 3 - 1) / + /</ 
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Figure 5. Left panel from I1UI : predicted nuclear modification for isoscalar- 
corrected 12 C, 56 Fe and 208 Pb to the neutrino- nucleus DIS stricture functions 
F^(xg,Q 2 ) (top) and xgF^(xg,Q 2 ) (bottom) versus Bjorken ij (left) and Q 2 
(right). The bands correspond to £ 2 = 0.09 - 0.12 GeV 2 . Right panel from fTH\ : 
Aqls calculated to 0(a s ) with charm mass (M c = 1.35 GeV) effects (stripes) 
and M c + resummed power corrections (band). Data is from CCFR, CHARM 
and IHEP-JINR. 



Similar rcscaling of Bjorken-x has been derived in for the EMC effect. 

• This resummation provides a natural framework to understand the differences 
in the nuclear shadowing of sea quarks, valence quarks and gluons. In 
the leading-order and leading twist parton model in v + A DIS F^(xb,Q 2 ) 
measures the valance quark number density with 4>vai.{ x ) ^ x~ Q " ai at small x. 
F^(xb,Q 2 ), a singlet distribution, is proportional to the momentum density 
of all interacting quark constituents and for xb <C 0.1 is dominated by the 
sea contribution, <j) sea {x) oc x~ aB ° a . Therefore, the xs-dependent shift from 
dynamical nuclear enhanced power corrections, Eq. I|ltjf> . generates different 
modification to F^(xb,Q 2 ) and F^(xb, Q 2 )- Predictions for the xb— and 
Q 2 — dependence of shadowing in v + A reactions is given in the left panel of 
Fig- El The right panel shows the high- twist contribution, to the Gross-Llewellyn 
Smith QCD sum rule 

S GLS - 3(1 - Aqls) = J dx B ^- (x B F^ A) + x B F^ A) ) . (17) 

• Gluon shadowing is also easily understood within the same resummation approach 
as the final state scattering of the struck gluon ^1] in p + A reactions. In the 
color singlet approximation gluons couple twice as strongly to the medium in 
comparison to quarks. The scale of power corrections then reads (C'a/Cf)^ 2 = 
2.25£ 2 = 0.20 — 0.27 GeV 2 and the corresponding dynamically generated gluon 
mass is twice as large. 

We now focus on the effects of dynamical nuclear shadowing on forward 
rapidity hadron production at RHIC. Any attenuation of the perturbative cross 
sections from coherent or inelastic parton scattering can be detected through 
the nuclear modification ratio R^ B , Eq. (fTTfl - In the special case of dihadron 
correlations such cross section reduction will be manifest in the attenuation of 




Figure 6. Left panel from 1141 : suppression of the single and double inclusive 
hadron production rates in d+Au reactions versus pt for rapidities y\ = 0, 1, 2, 3 
and 4. £ 2 = 0.12 GeV 2 . Also shown is the impact parameter dependence 
of the calculated nuclear modification for central 6 = 3 fm, minimum bias 
&min.bias = 5-6 fm and peripheral b = 6.9 fm collisions. The trigger hadron 
Pti = 1.5 GeV, yi = 3 and the associated hadron yi = 0. Right panel from 1141 : 
centrality dependence of 0*2 ( at moderate px, = 4 GeV, px 2 = 2 GeV and 
rapidities y\ =0,2 and y2 = 0. Central d+Au and p+p data from STAR. Also 
shown is C2(Ai/j) at small transverse momenta p^ = 1.5 GeV, Pt 2 = 1 GeV and 
rapidities y\ = 0, 4 and j/2 = 0. 



the area ^4f ar of the correlation function C2(A0). In order to compute the process 
dependent nuclear shadowing we isolate the small Xb dependence of single and 
double inclusive hadron production, Eqs. JSJ) and (JBJ), 

F a ^ cd (x b ) = \M ab ^ cd \ 2 . (18) 

x b 

The final state interactions of the struck parton lead to rescaling of the small 
momentum fraction Xf, as follows |14j : 

C2 



Fab^cd(xb) => Fabled Xb 



1 + C^A 1 / 3 - 1) 



(19) 



In Eq. (|19|l Ca keeps track of the representation of the struck parton with Cd = 1 
for quarks and Cd = 2.25 for gluon, as discussed above. 

Figure [S] show the predicted suppression of the single and double inclusive 
differential hadron production cross sections versus rapidity and centrality. The 
scale of power corrections was extracted for minimum bias collisions and the 
impact parameter dependence can be obtained via the corresponding rescaling 
with the nuclear thickness function. It is important to note that the value of 
£ 2 remains unchanged. However, in going from midrapidity to forward rapidity 
the perturbative scale —t that enters in Eq. 1|19[) changes from ~ 2pip to ~ pj,. 
This trivial kinematic dependence, which leads to larger suppression at forward 
rapidity, has been neglected in previous calculations. The right panel of Fig. 
shows the combined effect of the the small transverse momentum broadening 
and nuclear enhanced power corrections. In contrast to alternative hypotheses, 
the high-twist shadowing goes away quickly with the virtuality or transverse 
momentum as long as the large x a — > 1 threshold effects remain unimportant |36j . 
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Figure 7. Left panel from 1311 : suppression/enhancement ratio Raa{pt) for 
neutral pions at \fs NN = 17, 200, 5500 GeV. Solid (dashed) lines correspond to 
the smaller (larger) effective initial gluon rapidity densities at given y/s that drive 
parton energy loss. Right panel from |48|: predicted nuclear modification factor 
at y^ijvjv = 62 GeV for central Au + Au collisions. Enhancement, arising from 
transverse momentum diffusion in cold nuclear matter without final state energy 
loss is given for comparison. 



Further observations on the effects of nuclear shadowing on hadron production in 
d + Au at RHIC is available in |3J . Preliminary study of the effects of higher- 
twist resummation jHlEBEl on the low Q 2 modification to the DGLAP evolution 
equations has recently become available |38| . 

4. Jet quenching and jet tomography 

One of the few predicted |39j and observed signatures of nuclear dynamics associated 
with the presence of a hot and dense QCD plasma is the suppression of the high- 
er hadrons or jet quenching. While different attenuation mechanisms have been 
proposed HUJj it is the multiple inelastic parton scattering in the final state that 
was able to quantitatively predict many of the hadronic observables associated with 
jet quenching [221 1281 I42j . Different approximations to the medium-induced non- 
Abelian gluon bremsstrahlung dynamics can be found in |48| . In this section we give 
details of the Reaction Operator approach |44| . a momentum space iterative scheme 
for computing the radiation intensity. 

The full solution for the medium induced gluon radiation off jets produced in a 
hard collisions at early times r, et ~ 1 /E inside a nuclear medium of length L can be 
obtained to all orders in the correlations between the multiple scattering centers in 
the GLV approach 01]. The double differential bremsstrahlung intensity for gluons 
with momentum k = [xp + , k 2 /xp + , k] resulting from the sequential interactions of a 




Figure 8. Left panel from 1531 : the broadening of the away-side dihadron 
correlation function in central d+Au and Au+Au compared to scaled STAR data. 
In the bottom panel the broadening with and without suppression, <x Raa-, is 
shown. Right panel from |48l : attenuation of the double inclusive pion production 
cross section for pxi = Pt2 = 5 GeV at intermediate RHIC energies. 



fast parton with momentum p = [p + , 0, 0] can be written as: 



, dN g 
dx e? 2 k 



dN, 



in) 



n=l 



dx e? 2 k 



n=l 
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d?qi 



<5 2 (q t 



dAzi 



-2C 



(l,-,n) 



(m+l,---,n)(m,"-,n) 



COS 



COS 




(20) 



where ^\ = is understood. In the small angle eikonal limit x = k + /p + « uj/E. In 
Eq. i|2U|) the color current propagators, the formation times and the elastic scattering 
cross section are defined in 01] . For a derivation in the case of heavy quarks see 02] . 

Qualitatively, the behavior of the energy loss as a function of the density and the 
size of the system can be summarized to first order in opacity as follows: 

IE 



(AE) 



C R a s i? 
az — - — — z m 



2 A 



9C H 7ra~ 



pS(L) 2 In 



2£ 



9CRTTai 



dW 9 / r\ l„ 2£ 



■<L) In- 



static 
1 + ID 



(21) 



4 A ± dy \~/ JF(L) 

For static systems (AE) depends quadratically on the nuclear size. For the case of 
longitudinal Bjorken expansion this dependence is reduced to linear |47j but the energy 
loss is sensitive to the initial parton rapidity density dN 9 /dy. Understanding the 
effective color charge density dependence of AE the is the key to jet tomography 
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Figure 9. Left panel from 31: induced partonic multiplicity as a function of the 
experimental prcut for energetic quark jets at RHIC and the LHC. Right panel 
from 1541 : momentum density of hadrons associated with energetic back-to-back 
jets with and without medium induced bremsstrahlung. Secondary rescattering 
leads to gluon transverse momenta ~ 600 MeV. 



Transverse expansion does not affect significantly the attenuation of the inclusive 
spectra [17] . 

4-1- Applications of jet quenching 

The probability P(e) for fractional energy loss e = J^i^i/E due to multiple gluon 
emission is evaluated in the independent Poisson approximation [4fi| . One way of 
implementing the parton momentum attenuation is via the kinematic rescaling of the 
momentum fraction Zj — > Zj/(1 — e) in the decay functions, Eq. Q. Consequently 

D hj /j(zj,^) -> J <kP{e)^—D hj/j L^-,^ . (22) 

These corrections to the factorized formulas, Eqs. JSJ) and ©, can be easily 
implemented numerically. 

The left panel of Fig.|7|shows the predicted nuclear modification |3J as a function 
of ^snn at SPS, RHIC and the LHC as a function of the density of the quark- 

gluon plasma. The px dependence of R^a * s a resUi t of the interplay of the Cronin 
effect, jet quenching and nuclear shadowing |49| . The right panel shows a theoretical 
estimate for the 7r° (or tt + + tt~) attenuation at the intermediate RHIC energy of 
y?SNN = 62 GeV 0Hj • Such suppression is in agreement with the preliminary PHENIX 
and STAR measurements [231201 ior pr > 3 GeV. It is interesting to observe the better 
agreement between data and theory at the SPS with a recently extracted low energy 
p + p baseline |51| . Similar suppression at the intermediate RHIC energy has been 
found in |52|. 

For double inclusive hadron production qualitatively 1 < R^l A l R'aa 2 — 2- 
This attenuation will be manifest as a reduction of the area Ap al of C2(A</>), 
Eq. I|12|) . The left panel of Fig. [S] demonstrates that transverse momentum diffuion 
in Au + Au collisions is insufficient to reproduce the dissapearance of the away-side 
correlations [22] • Comparison to data require significant jet energy loss, compatible 
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with the attenuation in the single inclusives. The right panel shows the predicted 
quenching effect at ^/snn = 62 GeV 08] . 

The energy lost by the energetic partons during their propagation in dense nuclear 
matter is redistributed back in the partonic system. In the GLV approach 01] the 
medium induced virtuality is irradiated in higher frequency modes u> > oj p i ~ /i 
resulting in fewer more energetic gluons. The left panel of Fig. [21 shows the the 
recovered jet energy as a function of the experimental px cut and the corresponding 
induced parton multiplicity + l(jet). If the secondary (bremsstrahlung) gluons 
reinteract in the system [SI] their momentum is further degraded to pt — 600 MeV 
and 

1 l A£(r,A r) 
4 4 T(r,r) 

Numerical results for this scenario are shown in the right hand side of Fig. [5J For 
preliminary results on the recovery of the lost energy in jet measurements see |22II55| . 



N g (r,Ar) « 7 AS = - ; . (23) 



5. Conclusions 



The predictive power of perturbative QCD is based one the factorization approach that 
has been successfully tested in a variety of moderate and high transverse momentum 
processes |31|S1[7|E1E1E| m elementary N+N collisions. In this talk I discussed the 
theoretical and phcnomcnological aspects of the systematic incorporation of nuclear 
enhanced higher order and higher twist corrections in this formalism in £+A, p+A and 
A + A reactions. The current data on relativistic heavy ion reactions and the extended 
kinematic reach that is expected to become available due to improved statistics, RHIC 
upgrades, RHIC II and the LHC will help to further quantify the relative importance of 
the elastic j^HHSDl; inelastic 03] an d coherent [§1^11 El multiple parton scattering 
in cold and hot nuclear matter. 
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